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ABSTRACT 


To reduce the order of discrete-uncertain systems of higher order, a modified polynomial derivative approach is put forward 
in this study. In order to derive denominator and numerator of lower order model by means of affine arithmetic, a modified 
polynomial derivative stability array is built. To assess the respectivenumerator of the reduced order model and denominator of 
the reduced order model, Modified Polynomial Derivatives are used. The Calculation is based on the theories of Dolgin&Zehab 
[14]. In order to effectivelyassess the of the suggested strategy, it has been applied to two typical numerical examples. 
Additionally, it has been shown that the newer approach has a lower integral square error and offers a better approximation to 
the initial higher order system. 
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1. INTRODUCTION systems were presented by Bandyopadhyay et al. [12, 13, 


Low order models are often used to analyse and 
create complex systems, which are well acknowledged. 
The literature provides a number of approaches for 
scaling down large continuous time fixed systems. 
However, it is only possible to reduce discrete time fixed 
systems using a few numbers of approaches [1, 2, 3-11]. 

In the there 


approaches to reducing large-discrete time interval 


literature, haven't been many 
systems, which is true. Recent approaches involving 


model order reduction for large-discrete type interval 


14-16]. New concepts are proposed to improve the 
computation and output of the methods described below 
that may be found in the literature once they have been 
studied. 


2. PROBLEM FORMULATION 

These significant flaws and restrictions affect the 
approaches suggested by Bandyopadhyay and Ismail 
that are addressed [13]. Because this approach calls for 


the differentiation of the system's nth order polynomials, 
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it is computationally laborious. The approach sometimes 
generates unstable lower-order interval based models 
for stable original interval based discrete system of high 
order [13, 15], which is a serious drawback. To get 
denominator and numerator of the models with reduced 
order, a comparisonofthe initial 't' interval time instants 
of the reduced model and'm' Markov parameters from 
reduced model pared to those from the original 
uncertain model is done. This is done so that ttm = 2r. 
For this process, it is necessary to determine the 
time-instants of the original high order system G(z)and 
Markov parameters from the original high order system 
G(z), as well as those of the previously reduced rth 
order. This method has the ability to provide innovative, 
stable high order systems unstable models [13, 16]. With 
this method, the system's initial 't' interval time instants 
and'm' Markov parameters is matched with model that 
is resulting in t +m = r, and the denominator is obtained 
by maintaining the dominating poles of G. (z). The time 
basedinstants and Markov based parameters of the 
original discrete time system of higher order must be 
determined before using this approach. To determine the 
dominating poles in advance, the original system's 
denominator must be factorized [12]. The r-th order 
reduced model is created in a way that Rr(z) 
approximates G(z) to within 2r points using the Pade 
method. Transmutation that is reciprocal is necessary. 
Interval Routh type arrays must be created. It is 
necessary to identify the high order interval system's 
roots. A set involving 2r number of linear interval based 
equations must be resolved in order to produce a rth 
order reduced model [14]. 

The first "r" number of interval time instants from 
the model are to be matched with the first "r" interval 
time instants of the original system with high order to 
produce the numerator polynomial, and the dominant 
poles from the given discrete interval based system G(z) 
are used to generate the r-th order denominator Dr(z) of 
the reduced model. This approach is computationally 
expensive since the dominant poles must be determined 
from the denominator polynomial of the original 
system's denominator. The determination of the 
dominating poles is required for this procedure in 
advance. The essential goal and meaning of modelling 
are lost after a system's poles have been determined 
since its stability may then be evaluated directly. Similar 


to this, additional computational work is needed since 


the interval time based instants of the original system 
with higher order must be determined beforehand. 

Here we offer Modified Polynomial Derivative 
Technique, a novel and simple approach for order 
reduction of discretetime based interval systems of 
higher order, so as to overcome the drawbacks and limits 
discussed before. It is not dependent on long Routh type 
arrays, reciprocal transformations, Pade type arrays, or 
time-ahead calculations. The recommended method 
consistently yields "r-1" initial interval time moments in 
its reduced order model, improving the time response 
approximation, and low order interval models that are 
stable and unique. 

In this study, the previously created modified 
polynomial derivative method used for continuous and 
time based interval systems is expanded to the model 
reduction fordiscrete-time interval systemshigher order. 
The suggested technique is to be used for generating the 
reduced order models for the initial discrete time 
interval systems of higher order using the following 
procedural procedures. 


3. PROPOSED METHODOLOGY 
Give the initial high order discrete time interval 


based system the following form: 


To o Ez” +... 
CEIF = = m 
DE ED DO 2” Hsi ; 


A low order model is advised to be obtained for the 
original discrete time interval system of higher order 
discussed above. 

This is the definition for the suggested reduced 


order model with order "r" (r n): 


Reduction procedure: 
Step1:The original high order of transfer function 
from the p-domain will be derived as follows by 
applying a linear transformation, z = (p+1): 
G, (2) =G,(p)= 
z=(pt1) 


[BL Bp" +B} Be a] p" toe 


m m m-1? 


+[B, By] p+(By. Bo] 


C, (p) E i= + n = + n-l 
[A A ]p HIA, Apal Pt eee 


n-1? 


+[A A, ] p+[A).4g] 


ae 
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Step2: for FS V2 Bisa, Sab (r - 2) 


Using the suggested procedure, generate the model ESN (Bp) Be] (tera? ae 
; 3 [b; .b; 1=[B; .B; 1+ la; a; 1- [A; A; ] 
R:(p) of reduced order for aforementioned higher order [ADAG] 5 
system, Gn(p), as shown below. i-l fi At] [B7,B] | i 
Sse ed a Bo +5 CD ==], bl- - [a7 ,a7] 
ne n(p) [bb al p Pan +[b, .b, |p + [by bo | = [[A9 Ag | [By By] | 
(p)= = 
d(p) la, a} |p’ +la_,,a\,]p' + ies +[a; a; ] p+ [azsa Step5: 
--.(3) R:(w) is then transformed back into z-domain from the 


Step3:Reduced order denominator, di(p): linear transformation’s inverse, w = (z - 1), i.e. 


The following novel procedures are suggested to 


produce the reduced order model's low order £= R(w) 


w= (z-1) 


denominator polynomial aan By using the suggested strategy as in eqn, this yields the 


1 2: 
Forr=1, d,(p) = | ae Jian jee | Zi jazasi necessary lower order model (6). 


1 


Forr=2, 4. NUMERICAL EXAMPLES 
In order to assess the recommended- approach's 


; (im ey K, E 3k on adaptability and effectiveness, typical numerical 
per penal s eA: lee [40-40] instances are considered to apply to it, and the 


1 


3K 


3 


K 


1 1 


conclusions are successfully proven. 
Example 1: 


Consider the following stable discrete temporal system 
and in general, 


3 of fourth order [17]: 

ds n-r - r DOR K azr - j=1 
d (p)=| —— = iA; A lp + Sa aE vf 5 Eo 

"K rik 4 4 (@ = 

n-r j=l n-r D,(z) 
oe 4 2 
(4) [2.0,2.21] 2° - [4.424 4.911] z? + [3.3634 ,3.7334 ] z - [0.8721 ,0.968 ] 
2 P! P = 4 3 2 
where, Ko = ————— and el [1.0,1.05] z` = [3.037 ,3.189 ] z` + [3.425 ,3.596 ] z” - [1.6935 1.7782 ] z + [0.3084 ,0.3238 ] 
P!(P—@Q)! 


Step4: Redu chier Numerator, np): The following procedure is suggested for the above 


Na P heve been introduced for the original system in order to produce a 2nd order reduced 


derivation of the numerator polynomials with the first  nodel: 


interval Markov parameter and (r-1) first interval time Use of the suggested approach: 


based moments of the initial interval system. Stepl:Convert Giz) to Gs(p) using) a linear 


For ng, n = IIE) oly | transformation as shown below. 


For r = 2, n,(p) = [by 


„bi 1p +lby 561 


For r=3, ny(p) = ([b,.b;] p° + [b,.0, lp + [b-b] C (py) = a o) a Nn) 
z=(p+!1) D,(p) 


N,(p) 
ie, G (p)= — 
and in general, D, (p) 
3 2 
5 - r—1 = a 2 “y uana ~ 
en ae ee ae ee re, > ee BB [2.0,2.21](p +1) - [4.424 4.911 ](p +1)” + [3.3634 ,3.7334 |( p +1) 
- [0.8721 ,0.968 ] 
“ + [By By] = + = + = + = + 3 4 3 2 
where, [by bo 1= ———Lay ag] = [By By | Let (las aj] T [A43 I= 1; [1.0,1.05 ](p + 1)” - [3.037 3.189 ](p +1)° + [3.425 ,3.596 ](p +1)°] 
[Ay Ag ] 
une — [1.6935 1.7782 ](p + 1)+ [0.3084 ,0.3238 ] 
- + 
[b abala dla, a; ]; with d=mean of and 


RA [2.0,2.21] p? + [1.576 719] p’ + [0.5154 0.5414 ] p + [0.0644 ,0.0673 ] 


[1.0,1.05] p* + [0.963 1.011] p° + [0.314 ,0.329 ] p° + [0.0455 ,0.0468 ] p + [0.0026 ,0.0029 ] 


—Ė eee 
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Step2:In order to get a model R2(p) with lower order 


for aforementioned system, G(p), the following 


Step Response 


proposed SPDT is suggested: 


= 5 ‘oi th 
iss np) lb; bf] p +b; be] osl Original (4”order) J 
2 p = = - = 7 IL 

d,(p) lay va, le 40a, a; 194 bay a) = 


Reduced (2" order) 


Amplitude 


Reduced order denominator, d2(p): ost | 


: i 7 - 
d,(p)=[a},a}] p" +la,,a,]p+lay.ag] 


2 3 i 
K i : K l > r l , é 
=|- 2 114,,4,] p+ A 2 LAT, A ] p +A, Ag] 7 ; F; atsa ° 7 
K 2: K 2 

E" “Ks £y, K, 1 Fig. 1a: Step response of original 4torderand 274 
wit. = = a = ==? 
2 2 1 4 : 

tkaj 6 *K,} 2 order reducedsystem 


Thus 


d,(p) =[0.0523 ,0.0548 lp? + [0.0228 ,0.0234 ] p + [0.0026 ,0.0029 | 


Step Response 


Reduced order numerator, n(p): 


= ee EP Oe Original (4"order) 7] 
n,(p)=[b, ‚b; |p + [by .b9 1] 2 al | 
a Reduced (2™ order) 

0.3 - 4 

where , [b, „bi | =d parry |: oe | 

l [Bn Bn] [2.0,2.21] oa H | 
with d= mean of ————— = mean of ————— = 2.0574 | 
a H | 

[A, A, ] [1021705] o5 1 2 3 4 5 6 7 8 

Th Time (sec) 
us 


a _ 9 1 Fig. 1b: Step response of reduced 2”4 order and 
[b, „b; J= [0.1076 0.1127] and [by ,by ]= [By , Bo ]= [0.0644 ,0.0673 ] 
reduced 4 order systems 


Step3:Transform the reduced order model from the 
: r : X r Example 2: 
p-domain to the z-domain using the inverse linear 


transformation, as shown below: , , o. 
Consider the following stable 3rd order original 
n(z) 


R, (z)= R, (p) = 


= discrete time interval system: 
p=(@=1) d,(z) 


N (z) [3.25 3.35] z +[3.5,3.65] z + [2.8,3.0] 


[0.1076 ,0.1127 ] z — [0.0403 ,0.0483 ] G3(z)= = : 
R,(z)= s D(z) [5.4,5.5) 2° +[1.0,1.1] 27 +[1.5,1.6] z + [2.1,2.15] 
[0.0523 ,0.0548 ] z” — [0.0818 ,0.0862 ]z + [0.0324 ,0.034 ] 
(Stable)(Proposed Method) 


The step responses- from original discrete time 


The suggested technique yields the following 2nd 


interval system wh aaler and 2nd order model order reduced model for the aforementioned original 3"4 


. ; ._ order interval system: 
developed using the recommended technique is y 


[3.472 ,3.553 ] z + [6.078 ,6.447 ] 


compared in Figures 1a and 1b. R, (z) = - 
[5.733 5.867 ] z” + [1.667 ,1.799 ] z + [2.6,2.684 ] 


The recommended approach results in stable form 
(Using Affine Arithmetic) 


and reduced order models for a stable original discrete 


time interval based system of higher order, as shown in 
The 2”%4 order model obtained is: 
Figures 1a and 1b. 
[0.4717 ,0.5998 ] z + [0.9966 1.0075 ] 


Ry(z)= 
[1.4148 1.4589 Jz? — [0.8054 ,0.8674 ]z + [1.0,1.0] 


(Using Interval arithmetic) 
Step responses from the original discrete time 
interval system, the 2nd order reduced-models created 
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by recommended technique, and the interval arithmetic 
method are contrasted in Figures 2a and 2b. 


Lower bound 
2.5 T T T T T 


—— Original System 
AA 


= JA. 


Output 


o 1 1 1 L 1 1 1 1 1 
o 1 2 3 4 5 6 7 8 9 10 


Time (Secs) 


Fig. 2a: Step responses for original, AA and IA 
Upper bound 
— Original System 
— AA 
—— IA 


2.5 


Output 


L 1 1 
o 1 2 3 4 5 6 rå 8 g 10 
Time (Secs) 


i=} 


Fig. 
2b: Original Step reactions, AA, and IA responses are 


related to one another. 


5. CONCLUSION 
This 
derivative method for reducing the higher order discrete 


work presents a modified polynomial 
time interval based systems. Denominator polynomial 
and numerator polynomial are created from modified 
partial differentiation, respectively. Due to improved 
outcomes when compared to the current approaches, the 
suggested method's findings demonstrate its viability as 
a novel model order reduction method. The suggested 
Additionally, the 
suggested technique produces reduced levels of integral 


approach is readily calculable. 
square mistakes as compared to other already in use 


methods. 
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